arXiv: 1504.06318v3 [quant-ph] 28 Sep 2015 


Entanglement between exciton and mechanical modes via dissipation-induced coupling 

Eyob A. Sete^0 H. Eleuch^, and C.H. Raymond Ooi^ 

^Department of Electrical and Computer Engineering, 

University of California, Riverside, California 92521, USA 
^Department of Physics, McGill University, Montreal, Canada H3A 2T8 
^Department of Physics, University of Malaya, Kuala Lumpur 50603, Malaysia 
(Dated: September 30, 2015) 

We analyze the entanglement between two matter modes in a hybrid quantum system consisting 
of a microcavity, a quantum well, and a mechanical oscillator. Although the exciton mode in the 
quantum well and the mechanical oscillator are initially uncoupled, their interaction through the 
microcavity field results in an indirect exciton-mode-mechanical-mode coupling. We show that this 
coupling is a Fano-Agarwal-type coupling induced by the decay of the exciton and the mechanical 
modes caused by the leakage of photons through the microcavity to the environment. Using exper¬ 
imental parameters and for slowly varying microcavity field, we show that the generated coupling 
leads to an exciton-mode-mechanical-mode entanglement. The maximum entanglement is achieved 
at the avoided level crossing frequency, where the hybridization of the two modes is maximum. The 
entanglement is also robust against the phonon thermal bath temperature. 

PACS numbers: 42.50.Wk,03.65.Ud, 71.36.-l-c,78.67.De 


I. INTRODUCTION 

Hybrid quantum systems consisting of quantum me¬ 
chanical oscillators have become a platform for many in¬ 
teresting applications of quantum mechanics. In addi¬ 
tion to being a tool to understand the quantum to clas¬ 
sical transition, e.g., by creating entanglement between 
mechanical modes, mechanical oscillators have potential 
applications in quantum information processing. In this 
regard, there has been a growing effort in exploiting the 
mechanical degrees of freedom to engineer devices such 
as a microwave-to-optical (or vise versa) frequency con¬ 
verter and quantum memory ii. Moreover, quan¬ 
tum mechanical oscillator has been used as an interface 
to transfer a quantum state from an optical cavity to a 
microwave cavity [^, @ . Interest in merging optomechan¬ 
ical resonators with solid-state systems has been growing 
0; examples include, ultrastrong optomechanical cou¬ 
pling in GaAl vibrating disk resonator [sl, , cooling of 
phonons in a semiconductor membrane [lOl| , a strong op¬ 
tomechanical coupling in a vertical-cavity resonator |lll | , 
and surface-emitting laser [l2|. Coupling a mechanical 
oscillator to a microcavity that consists of a quantum well 
has been considered in the context of generating hybrid 
resonances [l^ among photons, excitons, and phonons, 
and studying the optical bistability [l^ [IJ . The physics 
of photon-exciton coupling has extensively been studied 
as presented in a recent review [T^. 

In this work, we analyze the entanglement between the 
mechanical mode and the exciton mode in a quantum well 
placed at the antinode of a microcavity that is formed by 
distritbuted Brag reflectors (DBRs). Even though the 
exciton and the mechanical modes are initially uncou- 
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pled, their interaction with a common quantized micro¬ 
cavity field results in an indirect coupling. We show that 
this coupling is a Fano-Agarwal-type [lilil coupling 
induced by the decay of the exciton and the mechani¬ 
cal mode caused by the leakage of photons through the 
microcavity to the environment (Purcell effect 1 1§. 1 1^ ). 
We analyze the entanglement in the adiabatic regime, 
where the damping rate of the microcavity exceeds the 
cavity-exciton coupling strength. A significant amount 
of entanglement between the exciton and the mechanical 
modes can be created at the exciton-mechanical mode 
hybrid resonance frequencies. We find that the maxi¬ 
mum entanglement between the two modes is achieved 
when the exciton and the mechanical modes hybridiza¬ 
tion is maximum. Surprisingly, the entanglement persists 
at high temperature of the phonon thermal bath. Our en¬ 
tanglement analysis is based on realistic parameters from 
a recent experiment [l^ . 


II. HAMILTONIAN AND LANGEVIN 
EQUATIONS 

We consider a microcavity formed by a set of dis¬ 
tributed Bragg reflector mirrors and consists of a quan¬ 
tum well placed at the antinode. The microcavity is cou¬ 
pled to the mechanical motion of the mirror via radiation 
pressure force and to the exciton mode in the quantum 
well. An exciton in the quantum well can be consid¬ 
ered as a quasi-particle resulting from the interaction be¬ 
tween one hole in the valence band and one electron in 
the conduction band. When the exciton radius is much 
smaller than the average distance between neighbouring 

_ 1 /Q 

excitons (^ riex ' with riex being exciton concentration), 
we treat the exciton as a composed boson. In general, 
in the weak excitation regime, where the density of the 
excitons is sufficiently low, the interaction between the 
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FIG. 1. Schematic of a microcavity made of two sets of dis¬ 
tributed Bragg reflectors (DBR) mirrors containing a quan¬ 
tum well (QW) and coupled to a mechanical motion (x) of 
the mirror. The quantum well is placed at the antinode of 
the microcavity so that the exciton-cavity mode coupling will 
be maximum. The black and the orange stripes correspond¬ 
ing to GaAs and AlAs layers, respectively. The microcavity is 
driven by a pump laser of normalized amplitude Sp and has a 
damping rate k. The DBRs are shifted from the equilibrium 
position due to the radiation pressure force. 


neighboring excitons due to Coulomb interaction is weak 
and can be neglected. However, in the moderated driv¬ 
ing regime, the interaction between neighbouring exci¬ 
tons becomes strong and nonlinear | 20 l - l^ . and leads to 
interesting properties such as squeezing and bistability 
(26l - l^ . In this paper we will consider the exciton as a 
composed boson. 

The coupled exciton-optomechanical system is de¬ 
scribed by the Hamiltonian 

H = LOaO^CL + ujexb^b + WmC^c -|- 

— goa^a(c + c^) -|- ig{a^b — ab^) + ab^b^bb. (1) 

Here the operators a, b, and c are annihilation opera¬ 
tors for a photon in the microcavity, an exciton in the 
quantum well, and a phonon in the mechanical oscilla¬ 
tor, respectively. The microcavity is driven by strong 
drive with frequency Wp; coa and Wex are respectively the 
bare microcavity and exciton frequencies. For the me¬ 
chanical oscillator, the resonance frequency is Wm and go 
is the single-photon optomechanical coupling; g is the lin¬ 
ear exciton-cavity mode coupling, and 2a = 6 e^aex/e^ 
[ 2 ^ is the nonlinear coefficient describing the exciton- 
exciton scattering due to Coulomb interaction with e, 
Oex) Cj and A being the electron charge, the exciton Bohr 
radius, the dielectric constant of the quantum well, and 
the quantization area, respectively. The strong drive of 
amplitude Ep = yjnP/huip with P and k being the drive 
laser power and the microcavity damping rate, respec¬ 
tively, leads to a large steady-state optical filed in the 
microcavity which increases the occupation numbers in 
each mode and the optomechanical coupling. The result¬ 
ing steady-state intracavity amplitude in turn shifts the 
equilibrium position of the mechanical oscillator through 
radiation pressure force. 

In the Hamiltonian (IT|) , the the first three terms in the 
first line represent the free energy of the system while the 


last term describes the coupling of laser drive with the mi¬ 
crocavity. In the second line, the first term describes the 
photon-phonon coupling, the second term represents the 
linear exciton-photon interaction, and the last term de¬ 
scribes the exciton-exciton scattering due to the Coulomb 
interaction. In a frame rotating with the drive frequency 
Wp, the interaction Hamiltonian m has the form 

V = — Aatt^a — Aexb^b + — goo^a[c -I- c^) 

-I- ig{a^b — ab^) + ab'^b^bb + iep{a^ — a), (2) 


where = LOp—oja, and Aex = ujp—uiex- Using the inter¬ 
action Hamiltonian Eq. we derive coupled equations 
for the macroscopic fields, a, 5, and c. These equations 
are obtained by replacing the operators with classical am¬ 
plitudes in the Heisenberg equations: 

a =——a + iAaO, + gb + igoa{c +c*) + Ep, (3) 
b =—^b + iAsJ> — ga — 2ia\b\‘^b, (4) 

c =-^c - iojmC + igo\af, (5) 


where 7 is the exciton spontaneous emission rate, and 
7 „i is the damping rate of the mechanical oscillator. The 
steady state solution to the above equations read 


^go\as\ 


Qs = -- 


“t” 

7/2 -|- i(Aex + 2 .alb) - 


bs = — 


gEp 


:7/4 -I- _ AoAex + f(KAex/2 -|- 7^a/2) ’ 


Aa{Ib)=Aa- 


( 7 / 2 )^ + A^^ 

+ (7m/2)2 g2 




Aex(db) — Aex + 2a/t,, 


( 6 ) 

( 7 ) 

( 8 ) 


where If, = | 6 sP is the steady state exciton number in 
the quantum well. Note that the Eq. (| 8 ]) yields nonlinear 
equation for in the form 


h 

9 ^ L 


fK-f 

V 4 


- _ A^Ae 


d' 


?A. 


— \^'n 


(9) 


The nonlinear equation for /{, is a signature that the ex¬ 
citon number can exhibit bistability [l^ . behaviour 
for a certain parameter regime. In the following, we 
discuss exciton-mechanical mode entanglement in the 
regime where the system is stable. 

The nonlinear quantum Langevin equations can be lin¬ 
earized by writing the operators as the sum of the steady 
state classical mean value plus a fluctuating quantum 
part: a = ttg -I- <5o, h = bs + Sb, and c = Cg + 6c. The lin¬ 
earized Langevin equations of the fluctuation operators 
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then read 

6a = ——5a + iAaSa + gSb + G{Sc + (5c^) + ^/nain, 

( 10 ) 

Sb = + iAex6b — g6a — 2iabl6b^ + V7^in, (11) 

6c = --^dc + iojm6c + G{6a^ — 6a) + •v/7mCini (12) 

where G = goV^ is the many-photon optomechanical 
coupling with hs = |dsp being the steady state mean 
photon number in the microcavity. For simplicity, we 
have chosen the phase of the coherent drive such that 
ds = —i|ds|. Here ain, hn, and c-m are the Langevin 
noise operators for the microcavity, exciton, and the me¬ 
chanical modes, respectively. All noise operators have 
zero mean, (ain(w)) = (6in(w)) = (cin(w)) = 0. We as¬ 
sume that the microcavity and the quantum well are cou¬ 
pled to a vacuum reservoir and thus the noise operator 
are delta-correlated: (ain(w)a;„(w')) = 2tt6{uj + uj') and 
(bin(uj)bl^(uj')) = 2 tt6{uj + uj'). However, the mechani¬ 
cal oscillator is coupled to a thermal bath and the noise 
operators have the following nonvanishing correlation 
properties in the frequency domain: (cin (uj)c\^(uj')) = 
27r(nth+l)^(w+w')> and (cl^(uj)cin(uj')) = 27rnth6(oj+uj'), 
where nth = [exp(/iuJm/kBT) — 1]“^ is the mean number 
of thermal phonons with ks being the Boltzmann con¬ 
stant and T the bath temperature. 


III. EXCITON-MODE-MECHANICAL-MODE 
ENTANGLEMENT 


We next study the entanglement between the exciton 
and the mechanical modes in the adiabatic regime, where 
the microcavity damping rate is larger than the exciton- 
cavity coupling, g, the cavity dynamics reaches qua¬ 
sistationary state. We then adiabatically eliminate the 
cavity mode degrees of freedom by setting dd = 0 in Eq. 
(HI. Substituting the resulting equation into Eqs. m 
and (HU, we obtain coupled equations for 6b and 6c that 
describe the dynamics of the exciton and the mechanical 
mode evolutions 


6b = - ^6b + i(Aex — 6uJex)6b — 2iab'^6b^ 

"^Gf)c{6c dc^) AfoCtin -j- 
6c = — — i{ujm + 6ujjn)6c — i5ujm5c^ — 

“f ~ ^cO-in + \/7mCin, 


(13) 

\Gbc6b 

(14) 


where F;, = 7 -I- 7 b with 7 b = Ag"^/ k[1 -I- (2Aa//c)^] be¬ 
ing the effective relaxation rate of the exciton due to 
the damping of photons through the microcavit y to the 
environment, also known as the Purcell effect |18l . 

Note that the relaxation rate of the exciton is increased 


by 7 b as result of interaction with the cavity mode; 
7 c = 4G^/«:[1 -|- (2Aa//c)^] is the effective damping rate 
of the mechanical mode. In contrast to the exciton 
mode evolution, the cavity-induced relation does not af¬ 
fect the decay term in the Langevin equation for 5c, it 
does however appear in the noise terms as manifested 
in Eq. (HI. Note also that the cavity-exciton coupling 
shifts the exciton frequency by dwex = 7 bAa/K. Simi¬ 
larly, the cavity-mechanical mode coupling gives rise to a 
shift dwni = 27 cAc/«: in the mechanical mode frequency; 
K{c) = y 7 K^(l-f 2 iAa/«:)[l-f ( 2 Aa/«;)^]“^/^ is the con¬ 
tribution of the cavity-induced dissipation to the noise 
operator of the exciton (mechanical) mode and finally 

Gbc = yjlblc (1 + 2iAa/Kj (15) 

is the effective exciton-mechanical mode cross coupling. 
Notice that the cross coupling depends on the effective 
decay rates 7 b and 7 c induced by the photon leakage 
through the microcavity, which is similar to the Eano- 
Agarwal effect [iElil. Dissipation-induced coupling has 
extensively been explored in quantum optics in creating 
coherence in three-level atomic systems |3ll - l33l| . Here we 
exploit the dissipation-induced coupling to entangle two 
matter modes: the exciton and the mechanical modes. 

To study the entanglement between the exciton and 
the mechanical modes, it is more convenient to use the 
quadrature operators defined by, 6xb = {6b'^ + 5b)/\/2, 
6yb = i{6b^ — 6h)l\f2, 6xc = (^c^ -I- 6c) 1^/2, and 5yc = 
i(^c^ — 6c) 1^/2 and similar definitions for fluctuation op¬ 
erators ajjqn, y_;,in (j = a,b). The equations for these 
quadrature operators in matrix form read 


u = Ru + g, 


(16) 


where u = {6xb,6yb,6xc,6yc)^ is vector of quadra¬ 
ture operators and g = with 

Fx,in — He(Ab)(rc^in F Im(Ab)ya,in T ■\/^^b,in^ Fy^^^ — 

R,e(Ab)ya,in Icn(Ab)(ra^in F \/^yb,\ni — ^7ni^c,in; 

and Fy^^^ = 2R,e(Ac)ya,in 2Im(Ac)(ra^in-t-^ 7 m 2 /c,in. The 
diffusion matrix R is given by 


R = 


A; 


\ -lm(Gbc) 


-A+x 

p + 

^ b 
2 
0 

—Re(Gbc) 


Re(Gbc) 0 ^ 

-lm(Gbc) 0 

~(Wm + 2(5(^ni) ) 


where F^ = Fb ± 4alm(6s) and A^ = Aex — jbAa/n F 
2aRe(6g). 

We focus on the steady-state entanglement between 
the exciton and the mechanical modes. For this, one 
needs to find a stable solution for Eq. dHl), so that it 
reaches a unique steady state independent of the ini¬ 
tial conditions. Since we have assumed ain, bm, and 
Cin to be zero-mean Gaussian noises and the correspond¬ 
ing equations for fluctuations bxj^m and 6yj^in are lin¬ 
earized, the quantnm steady state for fluctuations is sim¬ 
ply a zero-mean Gaussian state, which is fnlly charac¬ 
terized by a correlation matrix Vij = [{ui{co)uj{oo) F 






4 



FIG. 2. Logarithmic negativity En as a function of the 
detuning A = Aex = Aa normalized by the bare me¬ 
chanical resonance frequency cUra for the input laser power 
P = 24 fJ-W, and for various values of the thermal phonon 
number: nth = 70 (red solid curve), 100 (blue dashed curve), 
130 (green dotdashed curve), and 160 (magenta dotted curve). 
Here we use the experimental parameters from a recent ex¬ 
periment [T^ K = l/(5ps),7 = 1/(0.5ns),7m = l/(60ns), 
go/27r = 220 MHz,5/27r = 2.4 GHz, and tJm/27r = 20 GHz, 
and a = 10~^g. 

Mj(oo)Mi(oo))]/2. The solution to Eq. (fTCl) is stable and 
reaches the steady state when all of the eigenvalues of 
R have negative real parts. For all results presented in 
this work, the stability has been checked using the non¬ 
linear equation mentioned earlier. When the system is 
stable the correlation matrix satisfies Lyapunov equation 
RV + VR"^ = —D, where 

f ^ 0 0 0 \ 

]J — ^ ~2 ^ y/lblc 

0 0 ^(2rith + l) 0 

V 0 \/lblc 0 27c -b ^(2nth + 1) / 

and the elements of the drift matrix D are obtained using 
the correlations of the noise operators [sj defined earlier. 
Note that the cavity-induced dissipation terms contribute 
to the drift matrix. Notably, the off-diagonal element 
= Re(Gbc) contributes to the correlation between 
the exciton and the mechanical modes. 

In order to quantify the bipartite entanglement, we 
employ the logarithmic negativity E -^, a measure of bi¬ 
partite entanglement [351 136| . For continuous variables, 
En is defined as 

Fat = max[0,-ln2x], (17) 

where x = 2“^/^ [a — — 4detFj is the lowest sim¬ 

plistic eigenvalue of the partial transpose of the 4x4 cor¬ 
relation matrix V with a = det Va + det Vb — 2 det Vab 
[^ . Here Va and Vb represent the exciton and the 
mechanical modes, respectively, while Vab describes the 
correlation between the two modes. These matrices are 
elements of the 2x2 block form of the correlation matrix 

F = ( - I • The exciton and the mechanical 

V '^AB ) 



P(//W) 



P(//W) 

FIG. 3. (a) Logarithmic negativity as a function of the drive 
laser power P and for different values of the normalized de¬ 
tuning A/cJm = 1.05 (red solid curve), 1.10 (blue dashed 
curve), 1.15 (green dotdashed curve), and 1.20 (magenta 
dotted curve). Here we used the thermal photon number 
nth = 100. (b) Avoided level crossing between the eigenstates 
of the exciton-mechanical coupled system for rith = 100 and 
A/wm = 1.20. Notice that the maximum entanglement for 
A/wm = 1.20 in (a) appears at the power (P « 17.8 p W), 
where the maximum hybridization between the two modes 
occurs. All the other parameters are as in Fig. (2] 


modes are entangled when the logarithmic negativity En 
is positive. 

We numerically studied the exciton-mechanical mode 
entanglement by exploiting the indirect coupling medi¬ 
ated by the cavity field. Using realistic parameters from 
a recent microcavity experiment [l^ . we plot in Fig. [5] 
the logarithmic negativity Fat as a function of the nor¬ 
malized detuning A/uj-at and for different values of the 
thermal phonon occupation number, nth- Here we as¬ 
sumed the exciton-drive and microcavity-laser detuning 
are the same, = Aex = A. Figure [2] reveals that 
the exciton and the mechanical modes are strongly en¬ 
tangled, a demonstration of entanglement between two 
matter modes. The maximum entanglement is achieved 
at frequency where maximum hybridization between the 
two modes occurs. The entanglement expectedly de¬ 
creases when the thermal phonon number is increased; 
however, it persists up to thermal bath phonon number, 
nth < 200. 

In order to study the dependence of the generated en¬ 
tanglement on the applied input laser power, we plot in 
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FIG. 4. Logarithmic negativity as a function of the normal¬ 
ized detuning A/oJm optimized over the input laser power P 
range: 1-50 pW and for different values of the thermal phonon 
number: nth = 70 (red solid curve), 100 (blue dashed curve), 
and 130 (magenta dotted curve). 


Fig. [3] the logarithmic negativity versus power for dif¬ 
ferent values of the cavity-laser detuning. As can be 
seen from this figure, to obtain a maximum entangle¬ 
ment for a given cavity-laser detuning one has to apply 
a certain laser power strength. Naively, one would ex¬ 
pect that an increase in the coupling strength (due to 
an increase in power) to increase the entanglement. We 
however find that there exists an optimum amount of 
power that is needed to obtain the maximum entangle¬ 
ment for the realistic set of parameters [l2| • These peaks 
of the entanglement at different values of the laser power 
strength and detuning can be explained in terms of the 
exciton-mechanical mode hybrid resonances. The peaks 
appear at laser powers where the maximum repulsion be¬ 
tween the egienstates of the two modes occur [see, e.g.. 
Fig. 131(b)], indicating that the maximum entanglement 
is achieved at the maximum of hybridization. 

The optimized entanglement over the input power as a 
function of detuning and for different values of the ther¬ 
mal phonon number is shown in Fig. 01 The values of the 
cavity-laser detuning for which the peaks of the entan¬ 
glement occur shifts when the thermal phonon numbers 
are varied. This is because the effective coupling [see Eq. 
(USD] between the exciton and the mechanical mode de¬ 
pends on the cavity-induced damping rates. These damp¬ 
ing rates rely on the number of phonons, thus changing 
the resonance frequency at which maximum hybridiza¬ 
tion occurs. 

We note that the exciton-mechanical mode entangle¬ 
ment can be detected by measuring the optomechanical 
entanglement [38l - l40l| and the photon-exciton entangle¬ 
ment. From application view point, the generated entan¬ 
gled state has potential in one-way continuous-variable 


(CV) quantum computation [4l|. By forming a clus¬ 
ter of entangling gates, it is possible to implement CV 
quantum computation using our system. The exciton- 
mechanical mode entanglement might have advantages 
over that obtained between optical modes [4l| due to the 
robustness of the entanglement as well as the availability 
of semiconductor and micro-electromechanical (MEMS) 
technologies. The exciton-mechanical mode entangle¬ 
ment also means entanglement with mechanical oscillator 
or MEMS, a significant progress towards entanglement 
of macroscopic objects. Achieving entanglement in exci- 
tons against its large decoherence is an important step 
forward as it opens up new possibilities of merging quan¬ 
tum information with existing matured and ubiquitous 
technologies of semiconducting devices. 

IV. CONCLUSION 

In conclusion, we have analyzed the entanglement be¬ 
tween two matter modes (exciton and mechanical modes) 
in a hybrid quantum system consists of a microcavity, a 
quantum well, and a quantum mechanical oscillator. We 
have shown that although the exciton and the mechanical 
modes are initially uncoupled, their interaction with the 
common microcavity field results in dissipation-induced 
indirect coupling. This indirect coupling is responsible 
for the entanglement between the exciton and the me¬ 
chanical modes. Maximum entanglement is achieved in 
the adiabatic regime where the microcavity damping rate 
is larger than the coupling strengths and when the two 
modes form a complete hybridization. Recent successful 
experiments [1, [ifl, US] in coupling mechanical systems 
with microcavity pave the way for the realization of the 
proposed entanglement generation between exciton and 
the mechanical modes via dissipation-induced coupling. 
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